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Abstract 

It is well known that Raman scattering and fluorescence can be enhanced by the presence 
of metallic nanoparticles. Here, we derive simple equations to analyse the influence of metallic 
nanoparticles on upconversion processes such as non-radiative energy transfer or excited state ab- 
sorption. We compare the resulting expressions with the more familiar Raman and fluorescence 
cases, and flnd significant differences. We use numerical simulations to calculate the upconverted 
signal enhancement achievable by means of metallic spheres of different radii, and find particles 
of 100-400nm radius at infrared frequencies to be favorable. We also discuss the considerable 
challenges involved in using metallic particles to enhance upconversion for solar energy. 



I. INTRODUCTION 



Up- conversion processes allow to obtain photons of energy significantly larger than the 
excitation and are present, for example, in many different rare-earths systems. Potential 
applications include lasers [1], three dimensional imaging [2], quantum counters [3] or pho- 
tovoltaic cells [4, 5, 6, 7, 8]. An usual path to maximize the upconversion emission efficiency 
consists in trying different combinations of host matrices and rare-earth atoms. The crystal 
field from the host influences the transitions probabilities between the rare earth levels in a 
manner very dependent of the microscopic details of the quantum charge distributions. 

An alternative possibility, of interest here, relies on the possibility to locally enhance 
the excitation strength using metallic nanoparticles as local optical antennas. In addition, 
the decay rates of a given quantum emitter also depends on the environment, from purely 
electromagnetic, macroscopic considerations [9, 10, 11]. Metallic nanoparticles can thus be 
used as antennas both to enhance the incident field and to increase the emission rate. Surface 
Enhanced Raman Spectroscopy (SERS) serves as a clear demonstration of the potentiality 
of using such effects to obtain a strong signal. Beyond a chemical contribution, the strong 
local fields and radiative rates from plasmon resonances yield increases of many orders of 
magnitude. The first experimental demonstration of SERS was made with an ensemble of 
molecules adsorbed on a rough metal surface[12]. In the 1980s, as SERS experiments mainly 
concerned assemblies of emitters [13], no quantitative agreement with the theory could be 
found. The main reason is that the emission depends strongly on the orientation of the 
molecule, its distance to the antenna, and the exact shape of the latter. Recent advances on 
the manipulation at nanoscale enable the study of single molecule and single nanoparticles 
[14, 15]. SERS is now routinely used for sample characterization, and it is exploited in 
apertureless Scanning Near Field Microscopy to obtain images with subwavelength resolution 
[16, 17, 18]. 

A significant enhancement of the fluorescence from single molecules is also achievable[19, 
20, 21, 22, 23]. A recent review on the modification of single molecule fluorescence close to a 
nanostructure can be found in Ref . [24] . Nonetheless, the signal increase is in general much 
weaker than for SERS, due to a distinct light emission mechanism: Raman is a coherent 
process, comparable to Rayleigh scattering, but fluorescence is typically incoherent, with 
the population and depopulation of the different energy levels regulating photon emission. 
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In view of the fundamental differences between these emission mechanisms, it is not a 
priori obvious what will be the effect of metallic nanoparticles on upconversion processes. In 
2007, Polman et al. [25] showed that upconversion luminescence from Er ions could be used 
to image surface plasmons. The influence of small metallic nanoparticles on the upconverted 
signal from an assembly of Er ions has recently been investigated experimentally [26, 27], 
but no theoretical model was given. 

The purpose of this article is to address the issue of upconversion from a single emitter 
located in the near field of an individual nanoparticle. We consider two upconversion mech- 
anisms which are often predominant in systems consisting of rare earths such as Erbium 
or Erbium- Ytterbium: excited state absorption {ESA) and non-radiative energy transfer 
(NRET) [28]. One reason for optimism in the search of large signal enhancements is the 
proportionality between the upconverted signal far from saturation and the square (or a 
higher power) [29] of the light intensity. This reminds of the well-known fourth power rela- 
tion between the Raman signal and the local electric field enhancement. Thus, one could 
expect upconversion gains from localized plasmonic resonances as large as for SERS. Large 
enhancements are, however, far from guaranteed, as both NRET and ESA differ from Ra- 
man in an essential property. Like fluorescence, they are incoherent absorption-emission 
processes, which can be treated using rate equations, not a coherent scattering formalism. 
In other words, they involve real, not virtual, transitions. 

We derive simple analytical equations to study the influence of the environment on upcon- 
version processes. We also present similar equations for Raman scattering and fluorescence 
in order to compare the different processes. Our main objective is to obtain simple equations 
that captures the essence of the influence of nanoantennas on these processes. Thus, we use 
simplified level schemes and ignore some of the very complex physical phenomena that can 
enter in a real system. We illustrate the physical content of our model using numerical 
simulations of metallic particles of different size, at wavelengths where they exhibit localized 
Surface Plasmon Polariton (LSSP) resonances. We also briefly consider how differences on 
the level scheme translate into the enhancement of the upconversion signal (see Appendix 
B). 
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FIG. 1: (a) Scheme of the general scenario considered, where the emission of a quantum system 
is affected by closely situated particles. The incident and emission frequencies are different, (b-e) 
Energy levels used in most of this paper, for (b) Raman, (c) Fluorescence and upconversion via 
(d) excited state absorption and (e) non-radiative energy transfer 

II. MODEL AND NOTATION 

As Fig. 1(a) illustrates, we begin by studying a very general case, how particles imbedded 
in an homogeneous, loss-less medium affect the light emitted by quantum emitters with 
discrete energy levels - atoms, molecules, quantum dots... The energy level schemes in 
Fig. l(b-e) illustrate the different emission processes of interest, Raman, fluorescence, NRET 
and ESA. 

An experimental upconversion system can present a very complex energy level scheme, 
with the different transition rates -of possible electric dipolar, magnetic dipolar or even 
electric quadrupolar character- very sensitive to the exact microscopic details of the sample. 
To avoid considering so many variables, some of them not necessarily known, we consider 
only dipolar electric transitions and relatively simple schemes that illustrate many of the 
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phenomena of interest. We leave to Appendix B a further discussion of possible changes on 
the level scheme. 

We study the low intensity regime so that there are no saturation effects. We focus on 
the influence of purely electromagnetic effects that result directly from the presence of the 
metallic nanoparticles. We do not consider, for example, possible changes on the crystal 
field due to the presence of the particle. As we consider single, fixed quantum emitters, 
diffusion effects [28] which are often important for rare earths systems do not play a role. 
Non-local effects on the dielectric constant can in principle increase the non-radiative decay 
rate and thus diminish the emitted power. They are predominant at short distances to the 
substrate [30], but they will be significantly less important at the conditions of maximum 
upconversion enhancement shown in the following and we have not considered them. We 
also assume that the particle resonance is not modified by the presence of the rare earth 
emitters. 

For real systems, it will be necessary to match the particles resonances, the excitation 
frequency and the energy levels. To facilitate comparison over a broad range of excitation 
frequency ujin-, the energy difference between the ground and first excited state is always 
matched to the incident light, and other energy differences scale directly with this value. 
The excitation is assumed broad enough for population equations to be applicable and for 
frequency shifts [31, 32] due to the environment to be negligible, but sufficiently narrow to 
excite only the desired transition. The magnitude of interest is the ratio between the light 
power emitted in all directions at the desired frequency ojout with and without the particles 

Pout I ■ 

III. ANALYTICAL EXPRESSIONS 

We introduce in this section the equations describing the effect of the particule on the 
emitted power Pout/P^, for upconversion, both for ESA and NRET, as also for the known 
cases of fluorescence and Raman emission. Appendix A presents a more complete discussion 
on how to derive the equations and the corresponding assumptions. 
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A. Raman scattering 



We start with a short derivation of the well-known electromagnetic contribution to 
SERS[33]. Raman is a coherent process in which incident and scattered photons differ 
by a small energy amount that can be transferred, e.g. to a vibrational level of a molecule. 
The scattering enhancement is partly due to the enhancement of the local intensity at the 
scatterer position and partly due to the enhancement of the radiation emission due to the 
presence of the nanoantenna (appendix A). The former is characterized by the enhancement 
Kqi of the projection of the field into the dipole. The indices stands for the frequency 
huin = El ~ Eq of the field that matches a transition between two levels. The energy lev- 
els correspond to the excitation and emitted energy according to the simplified scheme in 
Fig. 1(b). The photon emission is characterized by the radiative rate Ffj in the presence 
of the nanoantenna and without antenna, at the emission frequency huout = Ei — Ef. 
The ratio of power emitted in all directions with and without nanoparticle is then given by 

P F^ 

- 'Aoil (1) 



We discuss in appendix A how in occasions l-ft'oil^ ~ ^fj/F™. The equality holds, for 
example, in the particular case of a confocal set-up. If the excitation and emission frequency 
are sufficiently close, we finally find 



out \ 1 4 



Kor\\ (2) 



po 

The often discussed proportionality between the intensity of the emitted light and the 
fourth power of the field enhancement follows from such frequently encountered conditions. 
It is a direct consequence of the coherent character of the process, not of the particularities 
of Raman, and also characterizes elastic (Rayleigh) scattering. 



B. Fluorescence 

Fluorescence is very different from the scattering the former involves a real tran- 

sition, and can thus be described by population equations. Introducing at this state the 
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notation for fluorescence and up-conversion, K^i gives as before the near field enhancement 
along the dipole orientation between levels |m> and |1>. F^j is the total decay rate between 
the levels and is the sum of three contributions F^j = (F^j + F^f - + /c^f ) that represent 
respectively the radiative decay resulting in photon emission, the non-radiative decay due 
to a transfer of energy to the particle and the non-radiative intrinsic decay due to internal 
transitions in the quantum emitter. Neither the process described by /c^f nor by F^^ result 
in photon emission to the far field, and the energy is lost, for example as heat. The quan- 
tum yield rj^^ = F^j/ ^mn) the fraction of quantum emitters in level |m> that decay 
directly to |1> by emitting a photon detectable in the far field at frequency uj^i. The sum 
over levels n serves to include formally cases in which the quantum emitter can decay from 
|m>to levels others than |1>. We note F^°, F^j and 7]^^ the values in the homogeneous 
medium, in the absence of any particle (Fj^f'°=0 by definition). 

In Fig. 1(c), a quantum emitter is excited to |1>, from which it decays to the fluorescent 
level |i> before returning to the ground state |0>. Appendix A shows that population equa- 
tions allow to account for the absorption and for the radiative and non-radiative processes 
and to derive the population of |i>. The quantum emitter decays from this level back to 
the ground state |0> with or without photon emission, as described by rj^Q. Finally, the 
comparison of the fluorescence signal with and without nanoparticle in the low intensity 
illumination regime under constant illumination yields the known equation [34, 35] 

Pout _ I 7>- 1 2 ^iO /Q\ 

~w - 

which only depends on the near-field enhancement and the ratio between the quantum yields 
with and without particles. A key difference with scattering is that for each excitation of 
the quantum emitter to |1> a maximum of one photon contribute to the fluorescence. Only 
after the quantum emitter has decayed a new excitation and subsequent photon emission 
can occur, and the number of the excited molecules is always lower than the total number 
of quantum emitters. Such limitation do not exist in scattering processes. If k^^^ = 0, the 
quantum yield for homogeneous, loss-less media rj'^Q is equal to one and the signal enhance- 
ment is lower than \Kqi\ , equal if the particles do not introduce non-radiative losses. 
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C. Up-conversion: Excited State Absorption 



Moving to discuss upconversion processes, we consider in this section the simple energy 
level schemes illustrated in Fig. l(d,e). In excited state absorption {ESA), a first photon 
excites a quantum emitter from level |0> to |1> and a second excites this quantum emitter 
to |2>. The energy of the emitted photon due to the |2> |0> transition is here exactly 
twice the corresponding value for the incident photons, which is not necessarily the case for 
more complicated level schemes. 

Proceeding similarly as before, we obtain in appendix A for the upconverted signal at 

^out 

J^out , 1 4 J- 10 '/20 

If = 0, the already introduced relationship ^p^s = \Koi\ is verified and there is no 

^ 10 

intrinsic losses, Pout/P^ < |-K^oi|^, equal in the absence of non-radiative decay. Under such 
conditions, the enhancement from ESA and fluorescence are similar. 



D. Up-conversion: Non- Radiative Energy Transfer 

The scheme in Fig. 1(e) represents the considered levels scheme for non-radiative energy 
transfer, NRET. The energy transfered non-radiatively from the relaxation of the two-level 
quantum system x according to the transiton |1> |0> serves to excite the three-level 
quantum system y from level |1> to |2>. 

Solving the equations in appendix A , we obtain for the upconverted signal from the 
transition |2> |0> of the quantum system y: 

- I |4 / ^ 10 \ ''20 /rX 

where we ignore any change in the efficiency of the non radiative energy transfer process 
due to the particle. We have also assumed identical field enhancement for x and y and 
an effective decay rate '^'^F^q for the complete system that is affected by the particle in the 
same way as for an individual quantum emitter. If there is no intrinsic losses, y decays 

2 ""^T^ 

directly from |2> to |0> and l-ft'oil = the presence of the particles can not increase 

^10 
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FIG. 2: A metallic nanoparticle influences the light emission from closely situated quantum 
systems, which can be oriented either radially or tangentially 

the emission of upconverted photons. 

We observe that the obtained expressions for the two upconverted processes are propor- 
tional to the fourth power of the near-field enhancement, similarly to the usual expression 
for Raman signal. However, the fourth power arises in upconversion because it is a two- 
photon process, while in Raman it emerges from a similar contribution from the near-field 
enhancement and the emission efficiency of the dipole in Eq. 1. Also, the upconversion 
processes involve real transitions as for fluorescence. Thus, when writing population equa- 
tions, non-radiative relaxation terms are introduced, and the interplay between intrinsic 
losses and radiative and non-radiative decay rates plays a key role. Due to these terms, we 
expect weaker enhancements for upconversion than for Raman scattering. In what follows, 
we shall analyse this interplay by modelling rigorously the electromagnetic interaction with 
a metallic nanoparticle. 

IV. RESULTS 

We study next light emission in the proximity of a spherical metallic particle embedded 
in glass (with dielectric constant egiass = 2.25). We consider (Fig. 2) dipolar excitation and 
emission parallel to the incident electric field vector and perpendicular to the propagation 
direction k with two orientations, one pointing toward the center of the particle (radial 
orientation) and the other tangential to the sphere (tangential orientation). The quantum 
emitter is in the plane parallel to and k and containing the center of the sphere, at a 
distance to the sphere surface d. For a sufficiently small particle, the total resulting electric 
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fields are exactly parallel to E'^ for both orientations of the dipole. For larger particles we 
still consider the transition to be only sensitive to the component parallel to 

We use numerical simulations to calculate the relevant parameters entering the power 
emission model: the enhancement factor and the modification of the radiative and non- 
radiative decay rates in presence of the particule. We then derive the modification of power 
radiated by the emitter in the different scenarios, Raman, Fluorescence and the two upcon- 
version processes. 

The simulations use a multiple multipole method (MAX-I) [36] to solve the Maxwell equa- 
tions. It is a semi-analytical method where the electromagnetic field is expanded by a series 
of basis fields called expansions. Each of the expansions is an analytical solution of the 
field equations within a homogeneous domain. The decomposition of the field is numerically 
optimized to minimize the error at the boundary conditions. The cross-section of the quan- 
tum emitter is considered small enough not to affect the near fields of the resonant sphere. 
The average error is smaller than 0.1 percent, often much smaller, and an ~ 25 percent 
decrease in the number of orders of the expansions does not significantly affect the results. 
The dielectric constants are taken from Palik [37]. 

A. Decay Rates and near-field enhancement 

Let us first focus on the influence of both the distance to the surface d and the wave- 
length on the parameters entering the analytical expressions for the emitted power: the 
field enhancement K and the radiative and non-radiative rates, F^/ F^° and r'^^/ F^°; the 
wavelength refers to the excitation, for the enhancement and to the emission, for the decay 
rates. 

Fig. 3 illustrates the distance and wavelength dependence of those parameters for the 
lOOnm radius silver sphere and radial orientation. It is useful for the discussion to also 
plot the quantum yield, with k'^^ = 0. The dipolar resonance is discernible in the near- 
field enhancement of ~ 4 at around ~ 1200nm(Fig. 3(a)). This wavelength is significantly 
shifted to the red of the small particle value due to the non-negligible size [38] . We can also 
observe at least one higher order resonance at shorter wavelengths (around 450 nm). The 
change on radiative decay rate F^/ F^° shown in Fig. 3(b) exhibits a similar behavior as 
the enhancement, with a broad maximum near 1200nm; in particular, near this maximum 
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FIG. 3: Dependence as a function of the wavelength and distance to the surface for a lOOnm 
radius silver tip and radial orientation of (a) the field enhancement (b) the change on radiative 
decay rate F^/ , (c) the change on non-radiative decay rate r^^/ and (d) the quantum 
yield for = 0. We note that the non-radiative decay rate (c) is plotted in logarithmic scale. 

pRy- pRo ^ 2 previously discussed. At shorter wavelengths the signature of not just 
one, but two higher order resonances are discernible. In contrast, the non-radiative decay 
rate change r^^/ F^*^ (Fig- 3 (c)) behaves very differently, with a very marked maximum 
near 350nm up to three orders of magnitude larger than the corresponding maximum for 
the radiative decay rate. F'^^/ F^'^ also increases significantly for small distances to the 
substrate. The conditions under which F^^/ F^*^ is significantly larger than F^/ F^° can 
clearly be seen on Fig. 3(d) when the quantum yield tends to zero. It happens for wavelengths 
around SSOram and at small distances. 

The size of the particle can greatly influence the results. This is illustrated in Fig. 4(a- 
c) where the wavelength dependence for changing radius and fixed distance to the surface 
d = 2nm are studied. If we concentrate on the dipolar resonance corresponding to the 
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FIG. 4: Dependence with wavelength of the (a,d) field enhancement K, (b,e) change on radiative 
decay rate pf'-/ r^-O g^j^^ ^g £-j ^.j-^g change on non-radiative decay rate pNRy'-pRO^ plots 
consider a distance to the substrate of 2nm. (a-c) correspond to silver spheres of different radii 
and radial orientation, (d-f) correspond to 100 nm radius, for silver and gold particles and radial 
orientation and for tangential orientation and silver particle 

first near-field enhancement and radiative decay rate maximum starting from the longer 
wavelengths, we observe a shift to the red for increasing radii. Furthermore, the peak 
becomes weaker and broader, with the maximum near-field enhancement remaining larger 
than 3.5 also for AOOnm radius particles. The non-radiative rate change behaves again 
differently than the near-field enhancement as seen in Fig. 4(c) , with the position of the 
maximum almost independent of size; it approximately corresponds to the real part of e^g 
equal to -e^; 

assi the resonance for a dipole in front of a —locally— fiat surface. 
The influence of the sphere material and orientation of the dipole are studied in Fig. 4(d-f) 
which represents the near-field enhancement and decay rates for a metallic particle of 100 nm 
radius of silver with radial and tangential orientation, and for a 100 nm gold particle for radial 
orientation. Whereas small gold spheres are known to present a dipolar resonance redshifted 
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FIG. 5: Dependence as a function of the excitation wavelength and distance to the surface for a 
lOOnm radius silver sphere and radial orientation of the change of emitted power Pout/ at the 
emission frequency of interest, (a) Raman signal, (b) fluorescence, (c) Excited state absorption and 
(d) Non-radiative energy transfer. The insert in (b) shows the fluorescent results for equivalent 
conditions but lOnm radius. 

with respect to silver, Fig. 4(d-f ) shows that the difference in position and strength is minimal 
for lOOnm radius; differences are nonetheless present for the observable high order resonances 
and the non-radiative decay rate peak. Finally, studying the tangential orientation for 
lOOnm silver particle leads to a completely different behavior of the enhancement and the 
radiative decay rate. Notably, the symmetry of the dipolar resonance implies strong electric 
fields and radiative decay rate for the radial, non tangential, orientation. 
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B. Signal enhancement 



Let us study now how the near-field enhancement and the modification of the radiative 
and non-radiative decay rates by the spherical particle affect the emitted power Pout/P^, 
for the different mechanisms of Raman scattering (Eq. 1), fluorescence (Eq. 3), ESA up- 
conversion process (Eq. 4), and NRET upconversion process (Eq. 5). To simplify, at this 
stage k^^ = for all involved transitions and, for ESA and NRET, the decay |2> |1> 
is negligible. To be able to study the evolution over a large wavelength range, we consider 
that the transition |0> |1> (and |1> —>■ |2>) coincides with the excitation frequency. 
The transition |1> |f> and |i> |0> corresponds to 0.95 and 0.8 times the excitation 
frequency for the Raman and the fluorescence level schemes, correspondingly. We discuss 
the results function of the excitation frequency. 

We start discussing the results for a lOOnm silver sphere and radial orientation (Fig. 5), if 
not otherwise mentioned. For Raman (Fig. 5 (a)), the presence of a simple metallic particle 
yields an enhancement of more than two orders of magnitude, with a broad maximum for the 
dipolar resonance and close distances to the particle. An enhancement of almost four orders 
of magnitude is possible at 2nm distance for lOnm radius. As expected, the Raman signal 
relates closely to the near-fleld enhancement (Figs. 3,4). An enhancement is also possible 
for fluorescence (Fig. 5 (b)), but signiflcantly weaker than for Raman, because it depends 
on the square of the near-fleld enhancement, not the fourth power, and because there is a 
competition between this enhancement and the quenching of the emitter. 

For the ESA case, as seen in Fig. 5(c), the signal is low at small distances to the surface d 
due to the quenching of the emitter. The signal is also specially low for excitation frequencies 
around ~ 350nm where the non-radiative decay rate for the |1> |0> predominates, and 
at around ~ 700nm, which coincides with strong non-radiative decay rate and thus low 
quantum yield at the upconverted frequency. Related to the discussion in the previous 
sections, we observe that the fluorescence and ESA enhancement look quite similar and 
both exhibit a maximum of 8-9 at the dipolar frequency. Considering a lOnm sphere (insert 
in Fig. 5(b) and Fig. 6(a)) serves to demonstrate that the similitude is not universal. 

The results are very different for NRET as seen in Fig. 5 (d). Indeed, the most striking 
feature for non-radiative energy transfer NRET is the very small maximum signal increase 
(only slightly larger than 1) in the presence of the particle. The increase appears at signiflcant 
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FIG. 6: Dependence as a function of the excitation wavelength and distance to the surface of the 
ESA signal Pout/P^- (a) lOnm radius, (b) 50nm radius, (c) 200nm and (d) 400nm radius silver 
spheres and radial orientation, (e) lOOnm radius gold sphere and radial orientation, (f) lOOnm 
radius silver sphere and tangential orientation. The insert in (f ) represents the results for identical 
conditions (included the wavelength and distance range of the plot) but for NRET. Notice that 
the figures (a-f) use different wavelength ranges. 

distances to the surface, i.e., where the particle influence is small. The dissimilarity with 
ESA signal is due to '^'^T\q/ ^'^T^q appearing squared for NRET. Indeed, this term describes 
how an increase of the decay rate of |1> can rapidly deplete this level and diminish the 
probability of an excitation to |2>, considerably affecting the upconverted signal. 

Fig. 6(a-d) and Fig. 5(c) illustrate the influence of the particle radius, for ESA and radial 
orientation. The main effect of small silver particles, of radius below 10 nm, is to decrease 
the upconverted signal, which can be attributed to strong non-radiative decay rates being 
predominant over near-field enhancement. For the larger particles considered (100 — 400r;,m 
radius), the value of the maximum upconversion enhancement is ~ 8, with the position of 
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the maximum shifting with the dipolar resonance, i.e. to the red with increasing radius. The 
optimal distance varies from ~ Ibnm to ~ bOnm for 100 and 400nm radius, respectively. 
The presence of higher order resonances and the influence of the non-radiative decay rates on 
r^o/ the |2> |0> quantum yield explain the presence of different secondary 

maxima for short wavelengths. For example, the clear minimum around 700nm for 50nm 
radius corresponds to an excitation frequency for which the non-radiative decay rate for the 
transition |2> |0> is particularly strong. 

Small spheres result in weak upconversion signal because at the position of strong en- 
hancement, the non-radiative decay rate can be very large. Further, the strongest 
near fields are obtained very close to the particle, where the quantum yield 7720 ^^^^ very 
small. For larger spheres, the larger spacial extent of the fields and the shift of the dipolar 
resonance means that sufficiently small non-radiative decay rates and strong enough near 
fields can be simultaneously obtained. We notice also that a signal enhancement has been 
achieved for small metallic particles[27, 39], but for a more complex scenario that makes 
understanding the underlying physics challenging. We discuss in Appendix B how more 
complicated level schemes than here can result in upconversion enhancement also for small 
particles. Further, centers such as pairs of silver atoms or ions are known to influence fluo- 
rescence via non-radiative energy transfers [40, 41] and could also play a role in upconversion 
experiments. 

Unsurprisingly from the results of Sec. IV A, the lOOnm gold spheres for radial orientation 
behave very similarly to their silver counterpart (Fig. 6(e)). Last, for the lOOnm silver 
sphere and tangential orientation (Fig. 6(f)), it is only possible to obtain a weak ESA 
signal enhancement of ~ 2. This reminds recent studies of fluorescence which point to the 
convenience of radial orientation [19]. In contrast, for NRET, wavelengths close to llSOnm 
and distances ~ 60nm we obtain an enhancement of ~ 3.5 for tangential orientation, larger 
than for radial orientation in the same conditions (Fig. 5(d) and insert in Fig. 6(f)). It 
is also larger than the values for ESA and tangential orientation (Fig. 6(f)), which is a 
consequence of the term F^q/F^q being larger than one, i.e., the total decay rate of the 
transition decreasing under those conditions. 

We have obtained a maximum enhancement of the upconverted signal for large particles 
on the ~ 100 — 400nm radius range, for wavelengths of ~ 1 to a few micrometers. The 
enhancement is about ~ 8 for ESA and radial orientation. The maximum position is at in- 
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frared frequencies and for a distance of a few tens of nanometers to the surface. Significantly, 
for large particles the maxima for r'^^/ and the near-field enhancement are not at the 
same frequency, which allows to favor one over the other. Changing the shape of the particle, 
or using several particles, should allow a similar effect also for smaller dimensions [42]. For 
NRET, we obtained a maximum enhancement of ~ 3.5 for tangential orientation. 

We emphasize that this study is mostly concerned by the emission of a single quantum 
emitter, which could be useful, for example, for near-field imaging or to characterize individ- 
ual hght sources. If the objective is to use upconversion for solar applications, the question 
of efficiency becomes crucial. In that case, one should also consider the fact that metallic 
particles not only enhance the field but also absorb the incident radiation. For example, 
we consider an absorption cross-section between levels |0> and |1> ctq^* = 10~^°cm~^, and 
that a 2 percent of the quantum emitters in level |1> are excited to |2> and result in an 
upconverted photon. In the presence of the resonant spheres and ignoring saturation, an 
enhancement of ~ 10 can be obtained resulting in an effective cross-section for upconversion 
~ 10 X 0.02 X 10~^° = 2 10"^^cm~^. This must be compared to the absorption cross-section of 
the particle. For a size parameter larger than one, we use the geometrical cross-section for a 
rough estimation of the scattering cross-section of a lOOnm radius sphere particle and obtain 
~ 7rl0~^°cm~^. The simulation results gives an absorption cross-section of ~ 2.5 lO^^^cm"^ 
at 1200nm. If we consider that all quantum emitters result in the maximum increase -no 
averaging over distance and orientation or considering interactions between the emitters and 
with the particles- it is still necessary to have 10^° rare earths atoms per particle to reach a 
similar number of upconverted photons than absorbed by the particle. For an homogenous 
distribution between and 100 nm, an unrealistic concentration of quantum emitters of 
~ 4 lO^^cm"^ is required. 

V. CONCLUSION 

To summarize, we have derived equations that give insight on how the emission from 
upconverted media behaves in the presence of nanoparticles, and we have compared them 
with fluorescence and Raman emission. Although the emission depends on the fourth power 
of the local field, the detailed analysis reveals typical gains much weaker than for Raman. 
It is sometimes similar but can also be significantly smaller than for fluorescence. We 
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obtain moderate enhancement of the upconverted signal when using large metallic particles 
illuminated in the infrared. We discuss the possibility of larger signal enhancements for 
more complex systems (see Appendix B). Using more complex resonant geometries can also 
be useful to obtain strong upconversion signal. For example, dielectric antennas can result 
in strong near fields without quenching [43]. 
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Appendix A : Derivation of the analytical expressions 
Raman scattering 

In this section, we derive the form of the enhancement of the scattered field produced 
by a nanoantenna. Let us consider a system where the excitation and the detection are 
produced by a source located at point ri. We characterize the source by an electric dipole 
moment pi = piu where u is a unit vector that describes the polarization of the source. 
The electric field produced at r reads : 



where G denotes the Green's tensor of the system at the excitation frequency ooin- This is 
simply the most general linear relation between two vectors. If the emitter is in a vacuum. 



then the Green's tensor is the vacuum Green's tensor Go- If the system includes a confocal 
microscope or any linear optical system, the Green's tensor accounts for it. Thus, although 
we have introduced a dipole source, a variety of illumination conditions is included formally. 
In the presence of a nanoantenna, the Green's tensor can be cast in the form: 




(6) 




(7) 
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where S accounts for the contribution of the nanoantenna to the Green's tensor. If we 
compare the field along v at r with and without the nanoantenna, we find an enhancement 
^u;i„ given by : 

v-G(r,ri,a),,Ju 

J^u>i„ = ^ (ej 

V • Go(r,ri,a;i„)u 

The field induces a dipole moment along v given by p2 = P2V = v(Q;eoV • E(r)), where a is 
the polarizability. It is either at the same frequency for Rayleigh scattering or at a slightly 
different frequency for Raman scattering cUout- The field produced by the induced dipole at 
a point r' is given by 

E(r') = G(r',r,cu„„t)p2 (9) 

We now consider two cases. The case where we have a confocal detection system and the 
case where we are interested in the emission over a given solid angle. Let us first assume 
that we use a confocal detection system with both illumination and detection at the same 
point r' = ri. We further assume that before detection light passes through a polarizer 
parallel to the unit vector u. The signal is then proportional to 



Pouto^ |u-E(ri)p = |u- G(ri,r,cj„„t)p2p 

P2 = Q;eov(v • G(r, ri, uJin)piu) (10) 

We now use the reciprocity theorem that implies that (^^^(r, ri, cUout) — Gnmi'^i, r, ujout)- It 
follows that Pout oc |v • G(r, ri, <^m)up|v • G(r, ri, <^out)up. If we now compare the signal 
Pout obtained in presence of the nanoantenna with the signal P^ without nanoantenna, we 
get : 

Pout ^ |v • G(r, ri, Uin)u\^\v ■ G(r, ri, o^onQup 
^ |v ■ Go(r, ri, cc;i„)u|2|v • Go(r, ri, UoutJW 
In other words, the enhancement of the emission is equivalent to the enhancement of the 
illumination. When the scattering process is elastic so that cuin — oUgut, we obtain an 
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enhancement given by 

|v ■ G(r,ri,cjj„)u| 



out 



|v ■ Go(r,ri,^^m)u|'' 



\K^J' (12) 



The equation also holds for Raman if the difference between emission and excitation fre- 
quency is sufficiently small. The exponent 4 shows that the signal varies like the fourth 
power of the enhancement produced by the presence of the antenna. The derivation shows 
that the fourth power enhancement is due to the role of the antenna that enhances both the 
illumination and the emission. Note that we have assumed that the detection is made at the 
same point and for the same polarization than the illumination. A different configuration 
might not yield this fourth power of the enhancement. 

In particular, we now consider that the detection system integrates over a solid angle 
n. We denote the emission rate of a dipole in the presence of the nanoantenna and 

^r-^° the emission rate without the nanoantenna. The emission rate refers here to the 
photons emitted through a surface subtending a solid angle Q. The enhancement of the 
Raman process is then given by 

p ^pR 

We have simply written here that the amplitude of the induced dipole is enhanced by a 
factor K and that its emission in a given solid angle is enhanced by a factor ^F^^^^. If we 
are interested in the emission over all directions, we simply use the corresponding radiative 
rate denoted by , which is the case in the text. Note that in the case that the emission 
pattern is not far from isotropic and is predominantly polarized as the excitation we have 

due to reciprocity \K\ ~ Q , a relationship that we have used in the text to obtain 

simplified equations; for ujin ~ ujout, the fourth power lif^.^l"^ is also under such conditions 
a good approximation for emission over all directions. 



Fluorescence scattering 

We present here the known derivation of the fluorescence equation[34, 35] to get familiar 
with the procedure followed for the upconversion expressions. As it involves a real transition, 
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the derivation used for the Raman case is not vahd. For weak illumination, the power emitted 
at the fluorescence frequency is due to spontaneous emission and is proportional to Ni, the 
population of the excited state. We write population equations for the different levels in 
order to derive the population of the excited state in stationary regime under constant 
illumination. The incident light excites the quantum emitter from state |0> to state |1>, 
with subsequent relaxation to |i> before decaying back to |0> (Fig. 1(c)). The different 
processes are absorption, radiative relaxation and non-radiative relaxation. As before, F^j 
describe the radiative decay rate between |m> and |1>, the non-radiative decay rate 
due to energy transfer to the particle and the decay rate from intrinsic losses of the 
quantum emitter [35]. The total decay rate F^j is the sum of the three. The absorption takes 
the form cr^fl^™' = 1/2 a'^^ ^/e/J^\K \ -E^i cos(6'i)p, with fi and e the permeability and 
permittivity of the medium, a^f the absorption cross-section for the concerned transition 
and the local intensity as deflned by the equality. cos(6'i) accounts for the projection of 
the fleld in the absence of particle to the direction of the dipolar transition of interest. 
A possible polarization rotation due to the particle is included in Kmi- The relaxation is 
considered to be fast enough for the population of the flrst excited level |1> to be negligible. 
Under the weak excitation assumption so that stimulated emission can be neglected, the 
rate equations are: 



N"^ refers to the population of level |m> and A^"^'' the population of all levels of the quantum 
emitter. A^"^'' ~ for weak excitation. The flrst term in the right hand side of the second 
equation describes the excitation of a quantum emitter from the ground level and the second 
term the total decay rate from the intermediate to the ground level at frequency Uio- 

In stationary regime, d/{dt) = 0. Since, Iq°{^"' oc lA'oip, we flnd from the population 
equations that A^' oc jA'oip/F^Q. The physics is clear : in the low intensity regime, the 
absorption increases linearly with the enhancement of the intensity l-K'oip so that the popu- 
lation of the excited state increases also linearly. Conversely, the increased total decay rate 
decreases the population. When writing that the emitted power is proportional to F^A,, 
we flnd that the power is proportional to the quantum yield rj^Q = T^/TJq. The comparison 



All 



dN° 
dt 




(14) 



21 



of fluorescence signals with and without nanoparticle finally yields the equation: 



Pout I 1 2 ^iO /-I p\ 

l^ = l^oi|^, (15) 
Up-conversion: Excited State Absorption 

Limiting ourselves for simplicity to the simple scheme in Fig. 1(d), the rate equations for 
ESA, under weak excitation intensity and neglecting direct excitation to level |2>, are 



All 



]\tO -.absf local 



01:^ + rTn + Tin 

01 



at . -10-- -20 



AtI --absjlocal 



dN'^ 0^12^12 at2tT Ar2 



- - r^o (16) 



Proceeding similarly as before, with A^^" ~ ^ A^-*- ^ A^^ due to the weak excita- 
tion, and V^^"^ = I'l"!"' for our energy levels and identical polarization dependence of both 
transitions, we obtain for the upconverted signal at 0^20: 

J^out I |4 10 '/20 

- ^01 1 -T^-Q- (I'J 



In general rj^Q = r2o/(l'2o + 1^21)5 but in the results section we considered to be neg- 
hgeable for both ESA and NRET. 

Up-conversion: Non- Radiative Energy Transfer 

Under the same assumptions than for ESA and the NRET scheme in Fig. 1(e), we can 
write 
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^ %1 _ -AE _ 

^ = ^iVi "iVi - ^N^ "Tl^ - ^N^ ^rfo (18) 

Where necessary to avoid confusion, a further index at the upper left side, as for example 
in '^A^-'-, distinguish between the quantum systems y and x, where y emits the photon at 
ujout and X is at the origin of the energy transfer. We have followed Page et al. [44] and 
assumed the two quantum systems in fast equilibrium, with an equilibrium decay rate ^'^F^q 
common for the the two levels |1> and the population relationship expressed by the fourth 
equality. There, is the partition function for the level m of x and AE (=0 in this 

paper) is the energy gap between the first excited level of y and x. The efficiency of the 
non-radiative transfer to level |2> is proportional to the population of both |1> levels, with 
0]^2 proportionality constant. We assume that the particle affects the equilibrium decay 
rate in the same way as for a simple quantum emitter. Solving the equations and using the 
definition of the quantum yield as before, we obtain for the upconverted signal 

^out i4 / J- 10 \ '/20 



where, assuming the two quantum systems sufficiently close, of same polarization dependence 
and with identical frequency difference between the ground and the first excited level, we 
have used l-ft'oil = K-K^oil = T-^oil- We have also neglected any change due to the particles 
in the energy transfer between the two quantum systems, otherwise a term of the form 
^12 would appear. 

Appendix B : Other Scenarios 

Real upconversion systems can have complicated level schemes, and their exact configura- 
tion may significantly affect the emitted signal. We illustrate now for ESA the consequences 
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FIG. 7: (a) ESA scheme, where the photon from |2> do not always decay directly to the ground 
level, (b-c) Excitation wavelength dependence for lOOnm radius silver sphere, radial orientation, 
15nm distance to the surface and the level scheme in (a), (b) Results for changing T^^j/ T^q, with 
Wont = l-9wi„, a;2iii = O.duJin and fc^f; = /c^g = k^^^ = 0. (c) Results for changing uj2iii/oJin, 
with LOout = 1.9iOin, rfP^i = rf\o and k^^, = kg§ = k^§ = . (d) Results for changing k^^/ TfP^, 
with a;2iii = 0.9u;j„, Uout = l-9u;j„, k^^§ = /c^f; = and r^j^; = 0. 

of several changes on the energy level scheme of the quantum emitter. 

Fig. 7 focuses on changes on the decay path for quantum emitters excited to level |2>. 
Such quantum emitters decay fast to an intermediate level |2i>, from which it can decay 
to |li> or |0> -emitting a photon, transferring the energy non-radiatively to the particle 
or losing its energy due to the intrinsic losses(Fig. 7(a)). The particular case |li> = |1> 
and |2i> = |2> reduces to the case illustrated in Fig. 1. For our numerical examples, if not 
otherwise mentioned, |li> = |1> and the |2i> |li> transition energy is 1.9 times the 
value for the excitation. The resulting equation looks very similar to Eq. 4 

Pout I 1 4 ^7210 /o^^ 

~ -'^oil -z^-Q— [20) 
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with the quantum yield of interest here taking the form 



^2i0 ~ pR I pNR 1 pR 1 pNR _i_ uNR _i_ l,NR (^■'■) 
^ 2i0 "T ^ 2i0 "T ^ 2ili ^ 2ili '^2i0 '^2ili 

Notice that intrinsic losses or the presence of the |2i> |li> always diminish the quantum 
yield term, but not necessarily the ratio of the quantum yield in the presence and absence 
of metallic particles [21]. Similar considerations apply to the decay 1/ F^q 

Fig. 7(b) concentrates on the influence of the strength of the decay path from level |2i> 
to |li>, with fc^fi = fc^^ = ky^Q = negligible for all transitions. At wavelengths around 
1400nm, slightly to the red of the dipolar resonance, the decay rate from the transition 
|2i> |li> is increased more strongly than the radiative decay rate of interest F^q, 
which negatively affects the upconverted signal. We notice that the effect is weaker for 
the wavelength of the maximum enhancement, because not only F2ii but also F^q increases 
significantly, the latter due to the effect of a higher order resonances. Indeed, for an adequate 
value of u!2m the increase of F^g can boost the quantum yield with respect to the value in 
the absence of particle; the enhancement of the upconverted signal is then larger than if the 
quantum emitter always decay directly to the ground state. The achieved improvement is 
moderate in our example(Fig. 7(c)). 

To analyse the effect of the intrinsic losses from the |2> level, we set F^^; = and 
change k^i^/ ^210- For particles where such intrinsic losses are considerable, an increase on 
the radiative decay rate F^q due the particle can lead to a considerable larger quantum 
yield and subsequent signal increase(Fig. 7(d)). 

Fig. 8 centers on changes on the decay path from the first excited state. We consider that 
a quantum emitter excited to level |1> decays fast and nonradiatively to |i> before either 
absorbing a second photon and being excited to |2> or decaying to the ground level. The 
resulting equation is again similar to Eq. 4 

p r™ n 

J- out I |4 ^ iO 720 /r,r,\ 

Intrinsic losses have a similar effect as before, and for large values the presence of the particle 
can significantly enhance the emission of upconverted signal (Fig. 8(d)) by increasing the 
term F^7 F^ = (F?^° + k^^^) / {V^ + F?^^ + A;/^^) . Using a Uio smaller than the excitation 
illumination allows to decouple the maximum of the total decay rate F^ and near-field 
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FIG. 8: ESA scheme, where the photon from |1> decay non-radiatively to an intermediate state 
before being excited by a second photon or decaying to the ground state, (b-d) Excitation avelength 
dependence for silver sphere, radial orientation, 15nm distance to the surface and the level scheme 
in (a), (b) Results for changing cjio/wm, with = k^Q^ = and lOOnm radius, (c) equivalent 
to (b) but for lOnm radius, (d) Results for changing k^^ / T^*^, for cuio = O.Owm, lOOnm radius 
and k^^ = 0. 

enhancement -ft'oij and thus to boost the achievable increase on the upconverted signal. 
The effect for lOOnm radius in Fig. 8(b) is relatively small because for the chosen conditions 
and the wavelengths around the maximum the total decay rate is not a strong function of uJiQ. 
We observe a much larger effect for lOnm radius, Fig. 8(c); non-negligible enhancement of 
the upconverted signal becomes then possible also for small particles, in contrast to Fig. 6(a) 
Copyright 2009 American Institute of Physics. This article may he downloaded for per- 
sonal use only. Any other use requires prior permission of the author and the American 
Institute of Physics. The following article appeared in Journal of Applied Physics and may 



26 



be found at R. Estehan et al., Journal of Applied Physics 105, 033107 (2009). 



[1] M.-F. Joubert, Opt. Mat. 11, 181 (1999). 

[2] E. Downing, L. Hesseling, J. Ralston, and R. Macfarlane, Science 273, 1185 (1996). 

[3] F. K. Fong, A. V. Lawrence, and R. E. de Wames, Phys. Rev. 170, 412 (1968). 

[4] A. Shalav, B. S. Richards, and T. Trupke, App. Phys. Lett. 86, 013505 (2005). 

[5] T. Trupke and M. A. Green, J. App. Phys 92, 4117 (2002). 

[6] G. Conibeer, Materials today 10, 42 (2007). 

[7] P. Gibart, F. Auzel, J. C. Guillaume, and K. Zahraman, Jpn. J. App. Phys. 35, 4401 (1996). 

[8] S. Ivanova, F. Pelle, A. Tkachuk, M.-F. Joubert, Y. Guyot, and V. P. Gapontzev, J. Lumin. 
128, 914 (2008). 

[9] H. Kuhn, J. Chem. Phys. 53, 101 (1970). 

[10] E. M. Purcell, Phys. Rev. 69, 681 (1946). 

[11] K. H. Drexhage, J. Lum. 1,2, 693 (1970). 

[12] M. Fleischmann, P. J. Hendra, and A. J. McQuillan, Chem. Phys. Lett. 26, 163 (1974). 

[13] A. Otto, I. Mrozek, H. Grabhorn, and W. Akemann, J.Phys.: Condens. Matt 4, 1143 (1992). 

[14] S. Nie and S. R. Emory, Science 275, 1102 (1997). 

[15] K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, I. Itzkan, R. R. Dasari, and M. S. Feld, 

Phys. Rev. Lett. 78, 1667 (1997). 

[16] M. Moskovits, Rev. Mod. Phys. 57, 783 (1985). 

[17] A. Campion and P. Kambhampati, Chem. Soc. Rev. 27, 241 (1998). 

[18] A. Hartschuh, E. J. Sanchez, X. S. Xie, and L. Novotny, Phys. Rev. Lett. 90, 095503 (2003). 

[19] M. Thomas, J.-J. Greffet, R. Carminati, and J. R. Arias-Gonzalez, App. Phys. Lett. 85, 3863 
(2004). 

[20] J. N. Farahani, D. Pohl, H.-J. Eisler, and B. Hecht, Phys. Rev. Lett. 95, 017402 (2005). 

[21] H. Mertens, A. F. Koenderink, and A. Polman, Phys. Rev. B 76, 115123 (2007). 

[22] P. Anger, P. Bharadwaj, and L. Novotny, Phys. Rev. Lett. 96, 113002 (2006). 

[23] S. Kiihn, U. Hakanson, L. Rogobete, and V. Sandoghdar, Phys. Rev. Lett. 97, 017402 (2006). 

[24] S. Kiihn, G. Mori, M. Agio, and V. Sandoghdar, Mol. Phys. 106, 893 (2008). 

[25] E. Verhagen, L. Kuipers, and A. Polman, Nano Lett. 7, 334 (2007). 

27 



[26] D. M. da Silva, L. R. P. Kassab, S. R. Liithi, C. B. de Araiijo, A. S. L. Gomes, and M. J. V. 

Bell, App. Phys. Lett. 90, 081913 (2007). 
[27] V. K. Rai, C. B. de Araiijo, Y. Ledemi, B. Bureau, M. Poulain, X. H. Zhang, and Y. Messad- 

deq, J. App. Phys 103, 103526 (2008). 
[28] F. Auzel, Chem. Rev. 104, 139 (2004). 

[29] M. Pollnau, D. R. Gamelin, S. R. Liithi, H. U. Giidel, and M. P. Hehlen, Phys. Rev. B 61, 
3337 (2000). 

[30] L A. Larkin, M. L Stockman, M. Achermann, and V. L Klimov, Phys. Rev. B 69, 121403 

(2004). 

[31] H. Morawitz, Phys. Rev. 187, 1792 (1969). 

[32] R. R. Chance, A. Prock, and R. Silbey, Phys. Rev. A. 12, 1448 (1975). 

[33] L. Novotny and B. Hecht, Principles of nano-optics (Cambridge University Press, Cambridge, 
2006). 

[34] J. Azoulay, A. Debarrc, A. Richard, and P. Tchenio, Europhys. Lett. 51, 374 (2000). 
[35] R. Carminati, J.-J. Greffet, C. Henkel, and J. M. Vigoureux, Opt. Commun. 261, 368 (2006). 
[36] C. Hafner, Post-Modem Electromagnetics: Using Intelligent Maxwell Solvers (Wiley, Chich- 
ester, 1999). 

[37] E. Palik, Handbook of optical constants of solids (Academic Press Inc., San Diego, 1985). 
[38] U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer, Berlin Heidelberg, 
1995). 

[39] V. K. Rai, L. d. S. Menezes, C. B. de Araiijo, L. R. P. Kassab, D. M. da Silva, and R. A. . 
Kobayashi, J. App. Phys 103, 093526 (2008). 

[40] C. Strohhofer and A. Polman, App. Phys. Lett. 81, 1414 (2002). 

[41] A. Martucci, M. de Nuntis, A. Ribaudo, M. Guglielmi, S. Padovani, F. Enrichi, G. Mattel, 

P. Mazzoldi, C. Sada, E. Trave, et al., App. Phys. A 80, 557 (2005). 
[42] P. Bharadwaj and L. Novotny, Opt. Exp. 15, 14266 (2007). 

[43] M. Laroche, S. Albaladejo, R. Carminati, and J. J. Saenz, Opt. Lett. 32, 2762 (2007). 
[44] R. H. Page, K. I. Schaffers, P. A. Waide, J. B. Tassano, S. A. Payne, W. F. Krupke, and 
W. K. Bischel, J. Opt. Soc. Am. B. 15, 996 (1998). 



28 



